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Abstract 

We consider the problem of triangulating a d-dimensional uniform grid of d-cubes 
that is cut by a ^-dimensional affine subspace. The goal is to obtain a triangulation 
with bounded aspect ratio. To achieve this goal, we allow some of the box faces near 
the affine subspace to be displaced. This problem has applications to finite element 
mesh generation. For general d and the bound on aspect ratio that we attain is 
double-exponential in d. For the important special case of d = 3, the aspect ratio 
bound is small enough that the technique is useful in practice. 


1 Introduction 

Recently the authors have developed the QMG mesh generator [10, 12], which is an algorithm 
for finite-element mesh generation of polyhedral regions in d dimensions. The main feature 
of the QMG mesh generator is that it guarantees an upper bound on the aspect ratio of all 
simplices occurring in the final mesh. The aspect ratio of a simplex is defined to be its 
maximum side-length divided by its minimum altitude. Small aspect ratio is important for 
accuracy in the finite element method [7]. 

The QMG algorithm generates boxes covering boundary faces of the polyhedral region 
P; these boxes are grouped into small clusters. Each such cluster contains boxes all of the 
same size and all the boxes in a cluster are associated with exactly one fc-dimensional affine 
subspace, which is a face of P. Thus, QMG gives rise to the subproblem of triangulating 
equally sized boxes cut by a hyperplane. 
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In this paper we consider this subproblem, namely, triangulating a uniform mesh of boxes 
in R d that is cut by a ^-dimensional affine subspace denoted F, where k lies between 0 and 
d — 1. We would like the triangulation to have the following properties: 

1. The triangulation is a valid simplicial complex, i.e., the intersection of any pair of 
simplices is a simplicial subface of both, and the union of the simplices is R d . 

2. The triangulation respects the affine subspace. In other words, the intersection of F 
with any simplex T in the triangulation must be a simplicial subface of T. 

3. The triangulation respects the mesh. That is, the intersection of any box B with T 
must be a simplicial subface of T. 

4. The triangulation has bounded aspect ratio. This means that there is an upper bound 
in terms of d and k only on the aspect ratio of all simplices. 

5. Every vertex of the triangulation is either a vertex of a box or lies on F. 

In general this problem may not have a solution because F may pass very close to a subface 
of a box. Therefore, we relax requirement 3 by allowing “warping” of boxes. We allow box 
faces close to F to be displaced. The amount of displacement must be bounded by a constant 
fraction of the box’s side length. 

For example, consider Figure 1 in which d = 2 and k = 1. Observe that the line F passes 
very close to the box vertices at (0, 1) and (1,2), and somewhat farther from (2, 3) and (3, 4), 
etc. Accordingly, we allow (0, 1), (1,2) and (2,3) to be displaced, and the boxes that contain 
them to be deformed in the course of triangulation. 

For the case d = 2 with k = 0, 1, warping rules were proposed by Bern, Eppstein 
and Gilbert [3] that solve this problem. In our earlier work [9] on mesh generation, we 
partially analyzed a warping rule for d = 3 using an argument based on cases. A case-based 
construction without an aspect ratio analysis was provided for d = 3 by Frey et al [6]. In 
this paper, we provide a complete analysis of the problem for all d and k. 

A different approach to guaranteed aspect ratio, which appears to be hard to generalize 
to dimensions higher than d = 2, was introduced by Chew [5]. 

2 Algorithm description 

To simplify notation, let us assume that the grid of boxes is the integer lattice; in other words, 
each box is of the form [si, Si + 1] X • • • X [sj, s^ + l] where «i, . . . , sj are all integers. A box face 
is a face of one of these boxes. A box face B can always be written as [si,ti] x • • • x [s^,^] 
where for each i, s, is an integer and t x is equal to either s, or s, + 1. We say that B extends 
over dimension i if the ith factor in this cross product is [s,, s, + 1] (rather than [s,-,s,]). The 
dimension of B, denoted dim(B), is the number of dimensions over which B extends. 
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Figure 1: In the d = 2, k = 1 case, the goal is to triangulate the mesh and the line crossing 
it while maintaining bounded aspect ratio. The result of our algorithm applied to this input 
is illustrated in Figure 2. 
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Our algorithm has two parts: preprocessing and triangulation. In the preprocessing 
phase, we identify the “close point” for each box face B. This close point of is a point 
lying on F that is close to B. 

We let cl(-B) denote the close point of B. Not all box faces have a close point; for a box 
face with no close point, we write cl(R) = 0. The rules for selecting a close point are as 
follows. 

First, assign a global priority order (a total order) to all box faces in the grid. Dimension- 
0 box faces (vertices) must have highest priority, dimension- 1 box faces next priority, and so 
on. Within a given dimension, the priority order can be arbitrary; for instance, lexicographic 
order on the coordinate entries is suitable. 

Next, we introduce parameters eo, . . . , td that specify how near F must pass to B in order 
to be considered “close” to B. We use a distinct parameter for each dimension, i.e., parameter 
e r is used whenever dim(£) = r. These parameters must satisfy several inequalities including 
the following: 

0 ^ Cd ^ ^ Ci ^ Co ^ 0.5. 

Further requirements on these parameters will be imposed below. In Section 7 we give the 
precise formula in terms of d and k for defining these parameters. 

Let B be a box face of dimension i, where 0 < i < d. Define N(B ) to be the infinity-norm 
neighborhood of size e,- around B\ in other words, 

N(B ) = {x : there exists y € B such that ||x — y||oo < e,}. 

Thus, for any B, N(B) is an axis-parallel parallelepiped. The rules for choosing the close 
point of B to F are as follows. We first assign close points to all dimension-0 faces, then 
dimension 1, and so on up through dimension d. The close point for a dimension-i face is 
determined by close points of its subfaces of dimension less than i if any exist. 

(Rl) If B has any proper subface with a close point, then among all such subfaces, choose 
the one B' with highest priority and let cl(£?) := cl (B'). 

(R2) If (Rl) fails, but the intersection of F with N(B) is nonempty, let cl(R) be the point 
in F fl N(B) closest to B. Here “closest” refers to infinity-norm distance. Break ties 
arbitrarily. 

(R3) If (Rl) and (R2) fail, then cl(R) := 0. 

Notice that closeness is inductive in the following sense. If B, B' are two box faces with 
B C B' and cl (B) ^ 0, then cl(B') ^ 0. Notice also that if the close point of B is determined 
by rule (Rl), and the particular subface with highest priority in rule (Rl) is denoted C, then 
the close subface of C must have been determined by rule (R2). The reason is that if cl(C) 
had been determined by (Rl) then C would have had a proper subface with a close point, and 
this proper subface would be a higher priority subface of B than C since lower dimensions 
have higher priority. 
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Now, with each box subface B we associate a point cl(B) defined as follows. If cl(B) ^ 0, 
we let cl(B) := cl(B). Else if cl(B) = 0, we let cl(B) be highest-priority dimension-0 subface 
of B. Thus, in all cases cl(B) is a point in R d on or near B. 

Now, finally, we explain the triangulation algorithm. Let a chain be a sequence of d + 1 
nested box faces Bo C B\ C • • • C Bj, where for each i, dim(B,) = i. Make a simplex 
that is the convex hull of cl(Bo), • • • , cl (Bd) and include this in the triangulation. Thus, the 
triangulation consists of all simplices induced by chains. The exception is the case of a chain 
for which cl(J3,) = cl(Bj) for some indices i ^ j. In this case the simplex associated with the 
chain is null and may be dropped. Notice that in the case of a box with no close subfaces, 
the rule in the last paragraph yields “Kuhn’s triangulation” [11] in which precisely d\ of the 
possible chains are non-null per box. 

For example, in Figure 2 we have illustrated the case d = 2,k = 1 with €o = 0.15 and 
= e 2 = 0. (The motivation for these values of eo, ei,e 2 is detailed in upcoming sections.) 
The infinity-norm neighborhood around each vertex is illustrated with a square. Notice that 
the line passes through the neighborhood of vertex (1,2). Thus, this vertex has a close point, 
as does every box edge and square that contains it. For instance, the square [0,1] x [1,2] 
has this vertex as a subface, and in fact, it is the highest-priority subface of this box with 
a close point. (The priority rule used in this figure favors objects with higher-numbered 
coordinates.) An example of a non-null chain in this box would be vertex (1,1) (with no 
close point), edge [0, 1] x {1} (with close point inherited from (0, 1)), and finally the square 
[0, 1] x [1,2] (with close point inherited from (1,2)). 

This concludes the description of our algorithm, except that we still need to select the 
parameters e 0 > • • • , Q. In the upcoming sections we describe the rules for these parameters 
and analyze the aspect ratio of our construction. 

3 Some of the e;’s are zero 

In this section we demonstrate two properties of (Rl) and (R2) that, without loss of general- 
ity, allow us to choose ed-k = £d-k+ i = • • • = = 0. Recall that k stands for the dimension 

of F and lies between 0 and d — 1. 

Lemma 1 Let B be a box face of dimension m such that m > d — k. Then rule (R2) can 
never apply to B. 

Remark. To show that rule (R2) never applies, we must prove that if F meets N(B), then 
F must also meet N(B') for some proper subface B' of B (so that (Rl) would take prece- 
dence). Consider Figure 2; consider specifically the intersection of F with a 2-dimensional 
box like [0, 1] x [1,2]. Notice that F cannot meet the neighborhood of this square unless it 
first penetrates the neighborhood of one of the vertices or edges of that square, which means 
that (Rl) would take precedence. This is the gist of the following proof. 

Proof. Suppose F meets N(B). Since B is a box face, we can write it as a product of 
intervals [si, <i] x • • • x [sd, td], where for each i, s, and t,- are integers. Furthermore, for each 
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Figure 2: The result of our algorithm applied to the geometry displayed in Figure 1. The 
dotted lines are the original grid. The neighborhoods around vertices are drawn as small 
squares. Each triangle of the triangulation has been slightly contracted in preparing this 
• figure. 
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i, either s, = U or si + 1 = t t ; let the set of indices in the latter category be denoted as J, 
that is, J is the set of indices over which B extends. Thus, | J\ = m. 

By definition, we know that 


N(B) — [^1 fmj d" ^m] X • • • X [sj Cm, tj T C m ]. 

Since F meets N(B), F is affine, and N(B) is a convex polytope, then their intersection 
F n N(B) is a nonempty convex polytope embedded in F. This polytope must have an 
extreme point, say x, and x must be the intersection of a face on the d — fc-dimensional 
skeleton of N(B ) with F. To simplify notation, let us suppose that this intersection point x 
lies on the following d — k face G of N(B): 

G — {si Cm} x {.S 2 c m } x ••• x {sfc c m } x [sfe+i e m , tfc+i d - c m ] X ••• x [sj e m , td -t- c m ]. 

Since m > d—k, there must be at least one index j € {1, . . . , k) that also lies in J. Consider 
the subface B' of B such that we replace [sj,tj] in the crossproduct of B by {sj}. This is 
a subface of B of dimension m — 1 containing G. Since e m _! > e m , we see that G defined 
above must be a subset of N(B'). Therefore, F meets N(B'). This shows that (Rl) would 
take precedence, so (R2) could not be applied to B. I 

The preceding lemma shows that, without loss of generality, we may take td-k+i = • • • = 

td = 0. 

Lemma 2 Let B be a box face of dimension exactly d—k and suppose that rule (R2) above 
applies to B. Then the close point of B in rule (R2) lies in B, i.e., its distance from B is 
zero. 

Remark. Consider the example in Figure 2. Consider the intersection of F with a neigh- 
borhood of a box edge like {3} x [3,4]. If F meets the neighborhood of this edge, but misses 
the neighborhoods of its two subfaces, then it must meet the edge itself. This is the gist of 
this lemma. 

Proof. To simplify notation, let us assume that B extends over the last d — k dimensions 
so that we can write it as 

B = {si} X ••• X {-Sfe} X [sfc+i , S)-+i d~ l] x ••• x [s^jS^d* l] . 

Let U be the affine hull of B , that is, the affine subspace of R d given by {si} x ••• x {s*} x 
R d-fc . Since F has dimension k , in the nondegenerate case F meets U at exactly one point 
x. We will assume nondegeneracy for now, and the degenerate case is considered below. If 
x G B then we are done. So suppose x £ B] we will derive a contradiction. 

We claim that x ^ B implies x £ N(B). To see this suppose x € N(B), i.e., for 
i = k + 1, . . . ,d, x,- € [s t - — td-k,Si + 1 + td-k}- (And we already know X{ = for i = 1, . . . , k 
since x 6 U .) Because x ^ B, there is at least one i 6 {k+ 1, . . . , d) such that x, ^ [s t , s, + 1], 
For instance, suppose x/t + i < so that Xk+i € [sfe+i — td-k, -sjt+i]. Then the bounds derived 
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on the entries of z in the three preceding sentences imply that x £ N(B') where B 1 is the 
subface of B given by 

B' = {sj} x ••• x { s k+1 } x [s k+2 ,s k +2 + 1] x ••• x [s^sj + l]. 

This is because e<i - k - 1 > td-k- But this contradicts the use of (R2) to find the close point of 
B. Thus, z £ N(B). 

Since (R2) applies to B, there is some point y £ F fl N(B). Consider the line segment 
from x to y. Since y £ N(B) and x N(B ), there must be a point z on the segment such 
that z £ N(B) but any point strictly between z and x is not in N(B). This point z lies 
on F because it is a convex combination of x and y. Furthermore, for each i = 1 
we know that z,- = s, and y,- £ [s, — td- k ,Si + td-k]- Therefore, every point on the segment 
from z to y, including z remains in the the interval [s, — td- k ,S{ + ej-fc] for all coordinate 
entries i = 1 Therefore, since z is a boundary crossing for N(B), there must be some 

coordinate i £ {1, . . . , k}, such that z t - is at the boundary of [s, — Cd- k , s, + 1 + ed - k ]: for 
instance, say that z k+i = s k+ i — td-k- But one checks that z £ N(B') where B' is the subface 
of B given above. Again, this contradicts the use of (R2). 

Thus, we see that the assumption z ^ B leads to a contradiction, so z £ B. 

We now consider the degenerate case, which also leads to a contradiction. Let y be the 
close point in F fl N(B ) determined by (R2). Because of the degenerate intersection of F 
and U, there exists a nonzero vector h £ R d such that the first k entries of h are zeros (i.e., 
h is parallel to U) and such that y + Xh £ F for all scalar values of A (i.e., h is parallel to F). 
Since N(B ) is bounded, for sufficiently large A, y + \h N(B). Thus, there is a choice of A 
such that y + \h £ N(B ) but y + \h ^ N(B) for all A > A. Since h has zeros in positions 
1, . . . , k, we know that y + \h must be at a boundary value of [s, — td- k , s, • + 1 + td-k] for 
some coordinate i not in {1, . . . , k}. But then, arguing as above, y + Xh would have to lie in 
N(B') for some subface B' of R, a contradiction to the use of (R2). I 

The preceding lemma shows that we can take td- k = 0 as well. 


4 Aspect ratio in terms of a determinant 


As mentioned in the introduction, the aspect ratio of a simplex is defined as its maximum 
side-length divided by its minimum altitude. Since all of the simplices in our triangulation 
lie inside a unit cube or a slightly distorted unit cube, the maximum side length is easily 
seen to be at most (1 + 2to)\/d. (In this section, distances are measured in the 2-norm.) 
Therefore, the aspect ratio is bounded above provided that we can show a lower bound on 
the minimum altitude. 

The purpose of this section is to obtain a lower bound on the minimum altitude in terms 
of a determinant. If S is a simplex whose vertices are vq . . . , Vd, let 


A = 



Vd - Vq 


( 1 ) 
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Although it appears that v 0 plays a special role in (1), in fact A is invariant under renum- 
bering of the i>j’s. Indeed, A is precisely d\ multiplied by the volume of S. 

Lemma 3 Let S be a d-simplex embedded in R d whose maximum side-length is s, and let 
A be defined by (1). Then the minimum altitude of S is at least A /s d ~ l . 

Proof. Let w be the altitude vector from Vd to the plane spanned by uo, . . . and 

without loss of generality, suppose it is the smallest. 

Then, since w is orthogonal to this plane, it satisfies the d — 1 equations { w , Vi — Vo) = 0, 
( w , v 2 — vo) = 0, . . . , ( w , Vd-i — vo) = 0, where (a, b) denotes the inner product of a and b. 
In addition, w satisfies the condition that its length is the same as the length from Vd to the 
plane; this condition can be written as (w,Vd — Vo) = (w,w). Let u be the vector that is 
parallel to w but satisfying {u,Vd — v 0 ) = 1. Thus, u is the solution to the square system of 
linear equations 


1 

CS 

1— l 

1 

cS 

O 


/°\ 

l>2 — ^0 


0 

• 

u — 

* 

0 
p 

1 

T-* 

1 

“O 


0 

L Vd-V o J 


l 1 / 


where u is regarded as a column vector. 

We would now like to obtain an upper bound on ||ti||. Let A denote the d x d coefficient 
matrix in (2). Note that (2) implies that u is the last column of A -1 . Therefore, ||u|| is equal 
to the (d, d) entry of A -r A -1 , that is, the inverse of AA T . The (d, d) entry of the inverse of 
AA T may be obtained by the adjoint formula for the inverse: 

|M| 2 = det(AA T )/ det(AA r ) (3) 

where A denotes the first d — 1 rows of A. We obtain an upper bound on the numerator 
by forming the product of its diagonal entries (recall that the determinant of a symmetric 
positive definite matrix is bounded above by the product of its diagonal entries). The 
ith diagonal entry of AA T is the square of the length of the edge from i \ to uo- Thus, 
det(AA T ) < s 2 ^ -1 ). Next, we need a lower bound on the denominator of (3), which is equal 
to det(A) 2 . We see that A = | det( A) | by (1). Thus, the denominator of (3) is A 2 . 

Combining all this, we conclude that ||u|| < s d-1 /A. Recall that w is a multiple of u 
scaled so that ( w , Vd — v 0 ) = ( w , w). Because u satisfies ( u , Vd — u 0 ) = 1, it is easily seen that 
we must take w = u/||tt|| 2 . Thus, ||tt>|| = l/||u|j, hence ||iw|| > A/s'* -1 . I 

5 Bounds on the coordinate entries of close points 

In this section we begin the analysis of the aspect ratio. The purpose of this section is 
to obtain bounds on the coordinate entries of the vertices that occur in a non-null chain. 
These bounds are used in later sections to obtain constraints that the e;’s must satisfy. The 
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Figure 3: This figure shows what goes wrong if eo is too small. In this figure, d = 2 and 
k = 1, the edges of the square are dashed, and the eo-neighborhoods of the vertices are the 
small squares. If e 0 is very small as in the figure, then we get a triangle with bad aspect 
ratio, namely, the skinny triangle chopped off the left edge of the square. Thus, we want e, 
for each i = 0, . . . , d — k — 1 to be as large as possible subject to other constraints. 
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(0,1) 


(1.1) 



Figure 4: This figure shows what goes wrong if eo is too large. In the above example of the 
d = 2, k = 1 case, if to is too large, then we can get a chain with three different close points 
(marked with asterisks) lying on F. For instance, in the above figure we see that F has a 
close point for vertex (0,0), a close point for vertex (1,0) which is also the close point for 
edge [0, 1] x {0} (assuming a certain priority rule), and finally a close point for vertex (0, 1) 
which is the close point for the whole square. Thus, we would try to form a simplex that is 
the convex hull of these three distinct close points. This simplex would have infinite aspect 
ratio because the points are colinear. Notice that this simplex is not “null”; it does not have 
a repeated vertex. Therefore, it cannot be dropped from the triangulation because the result 
would no longer be a simplicial complex. It is clear from the figure that we can prevent this 
difficulty by making e 0 sufficiently small with respect 1 — e 0 . 
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Figure 5: This figure shows what goes wrong if is too large with respect to e 0 in the 
d = 3, k = 1 case. The figure shows a planar projection of the 3-cube. The neighborhoods 
around vertices and edges are drawn with t\ = eo- The set F, a line, is perpendicular 
to the plane of the page at the point marked with an asterisk. Notice that F meets the 
neighborhoods of three edges of cube marked a, b, c but misses the neighborhood of every 
vertex. This results in a chain whose associated simplex has infinite aspect ratio as in 
Figure 4. In order to prevent this, we must make to sufficiently larger than ei. 
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reader may wish to consult Figure 3, Figure 4 and Figure 5 for informal examples of these 
constraints before reading these sections. 

The main idea of the bounds in this section can be summarized as follows. Recall that a 
chain is a sequence of nested box faces. Therefore, a face B x in the chain extends over one 
additional dimension, say dimension j, compared to In this section we show that the 

close point for B, must have a jth coordinate entry that is a sizeable distance from £,•_ j, 
else the simplex would be null. 

Let us introduce notation for this section and the upcoming sections. Suppose that 
Bo C Bi C • • • C Bd form a chain of box faces that determine non- null simplex, i.e., a 
simplex with no repeated vertices. By reflecting, translating, and permuting indices, we can 
assume without loss of generality (to simplify notation) that the chain is 

B 0 = {0} x {0} x ••• x {0}, 

Bi = [0,1] x {0} x ••• x {0}, 

B 2 = [0, 1] x [0, 1] x {0} x • • • x {0}, 

B d = [0,1] x ••• x [0,1]. 

For i = 0, . . . ,d, let i>, be cl (£?,•). The assumption that the simplex is non-null means 
that the u,- are all distinct. Let / be the lowest index such that Bi has a close point. Then 
all of Bi, • • • , B d have a close point by the inductiveness of (Rl). (Let / = d + 1 if no box in 
the chain has a close point.) 

Assume the u,-’s are written as row vectors, and let V be the (d + 1) x d matrix 


v 0 


V = 


v d 


Let the rows of V be numbered 0 to d, and the columns 1 to d. The purpose of this section 
is to derive some bounds on the entries of V. 

We start by looking at the top l rows of V. Note that each of these rows is a box vertex, 
hence all the entries in these rows are 0’s and l’s. Consider row i of V for i = 0, 1 
(assuming / > 0). Row i is a box vertex lying on Bi. For any j = i + l,...,d, Vij = 0 
because all points in Bi have zeros in positions i + 1 , . . . , d. 

We claim further that for i = 1, — 1, Va = 1. Suppose not; suppose Vu = 0. This 
means that n, is a vertex in R,_i as well as in Bi. But recall that v, is the highest priority 
vertex in Bi. Since B t is a superset of Bi- 1 , we conclude that v, must also be the highest 
priority vertex in Bi-\. This in turn implies that v, = n,_i, contrary to our assumption that 
the simplex is non-null. This proves that Va = 1. This concludes the analysis of the first l 
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rows of V . For example, if / = d + 1 then V has the form 


' 0 0 O' 

1 0 0 

x 1 0 0 

x x 1 0 

x x 1 


where each ‘x’ is either 0 or 1. 

Now we consider rows of V numbered l to d. Recall that each of these rows has a close 
point. This close point was determined by either (Rl) or (R2) in each case. Let Ci, . . . ,Cd 
be subfaces of Bi,...,Bd defined as follows. If (Rl) was applied to Bi, then let C, be the 
high-priority proper subface of B, that determines v,\ If (R2) was applied to B{, then let 
Ci := Bi. Thus, in either case, u, is the close point of (7,, and Vi € N(Ci). 

Let the dimension of Ci, for i = /, . . . ,d be denoted as r(i). Note that r(i) < d — k for 
any i as proved in Section 3. We claim that r(l) > r(l + 1) > • • • > r(d). The reason is as 
follows. Consider an i € d — 1}. Notice that Ci, as a subface of Bi, is also a subface 

of Bi+i. Therefore, C t+ i must have higher or equal priority compared to C, since C, + i was 
chosen as the high-priority subface of R,+i. Thus, C,+i must have lower or equal dimension 
than Ci because lower dimensional faces have higher priority. 

Now, fix a particular i between / and d. We claim that \Vij\ < e r (q for each j = i+1, . . . , d. 
This because v, lies in N(Ci), and C, is a subface of Bi. Therefore, all points of Ci are zero 
in positions i + 1,. .. ,d. Thus, a point in N(Ci ) has absolute value at most e r (,j in these 
positions. 

Now consider the diagonal entry Vu for some i > 0 between l and d. There are two cases; 
either r(i ) = 0 or r(i ) > 0. If r(i) = 0 then Ci is a box vertex. By the same argument as 
above, this vertex must have a ‘1’ in the ith position (otherwise u, would also be the close 
point for B,_i). Therefore, the ith entry of u, is at least 1 — cq. 

The other case is r(i ) > 0. In this case, we claim that Vi,- > e r (,-)_x. There are three 
subcases: in the ith coordinate, either C,- is identically 0, identically 1, or extends over 
dimension i. 

If Ci is identically 0 in the ith coordinate, then C, is actually a subface of R,_i . But 
then, as above, it must be the highest priority subface of R,_i, contradicting the assumption 
that the simplex is not null. So this case is impossible. 

If Ci is identically 1 in the ith coordinate, then every point x 6 N(Ci) has x,- > 1 — e r (,) > 
0.5 because all of the e,’s are less than 0.5. Thus, V^, > 0.5 > £ r (t')-i- 

If Ci extends over the ith coordinate, then it has a proper subface C- obtained by fixing 
its ith coordinate to zero. The dimension of (7,- is r(i ) — 1. Notice that u, cannot lie in N(C'i ) 
because (7 t - has higher priority than C, and is also a subface of Bi. But notice that if V t , 
were e r (.)-i or smaller, then indeed u, would lie in N(C^) because e r (i)-i > e r (,). Therefore, 
to avoid this contradiction, we must have Vu > e r (i)-i- 
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To unify this case with the case of r(i) = 0, let us adopt the notational convention that 
e_i = 1 — e 0 . Then we can say that in all cases, V;,- > e r (*)— i for * = max(/, 1), . . . , d. 

Finally, we would like to estimate entries of V below the diagonal. The crude upper 
bound that 

|^«j| 5- 1 + Co 

follows because all the v^s lie within e 0 of the unit cube. 

This concludes our estimates of the entries of V. To summarize: 

• The “superdiagonal” entries, that is, entries Vij with j > i, satisfy: 

- Vij = 0 if i < /; 

- | Vij I < e r(0 if i > l. 

• The “diagonal” entries, that is, entries Vu for i > 1, satisfy: 

- Vu = 1 if * < /; 

- Vu > e r(t')-l if i>l- 

• The “subdiagonal” entries, that is, entries Vij with j < i, satisfy | | < 1 + eo. 

Finally, we have showed that r(l) > r(l + 1) > • • • > r(d). 


6 A bound on the number of close points 

In this section we show that l > d — Jc, i.e., at most k + 1 of the boxes in the chain can have 
close points. To prove this requires bounds on the e,’s which are derived below. (Recall the 
example of Figure 4, in which d = 2, k = 1, and 1 = 0. The result of this section is that l > 1 
must hold for the d = 2, k = 1 case provided Co is chosen correctly.) 

Lemma 4 Suppose the sequence of e,- ’s satisfies (12) below. Then at most k + 1 box faces in 
the chain can have close points (assuming, as in Section 5, that no close points in the chain 
are repeated). 


Proof. Suppose not; suppose that the last k + 2 boxes, that is, boxes Bd-k-i,' • • all 
have close points Vd-k- 1 , . . . ,Vd- Consider forming the {k + 1) x d matrix 


Vd-k — Vd-k - 1 
Vd-k+1 — v d-k - 1 


Vd ~ Vd-k - 1 


Notice that each row of Y is the difference between two points lying on F; therefore each 
row lies in a ^-dimensional vector space. Thus, the maximum possible rank of Y is k. Define 
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Y to be the submatrix formed by the last k + 1 columns of Y . Thus, Y is a (k + 1) x (k + 1) 
matrix whose rank is also at most k. This means that the determinant of Y is zero. We will 
show that for appropriate choices of the e,’s, it is impossible for this determinant to be zero, 
yielding a contradiction. 

Let us number both the rows and columns of Y with the index sequence d — k, . . . ,d 
since this is the natural index sequence inherited from V. Recall that the determinant of Y 
is a sum of ( k + 1)! terms. We show that the term from the diagonal must be positive and 
sufficiently large that the combination of all the other terms could not cancel it out. 

Let this diagonal term be denoted as 9: 

9 = Y d - k , d - k ■ ■ ■ Y d<i . (4) 

Recall that Yu = Va — V d -k-i,i • Therefore, using the bounds from the last section, we 
conclude that 

Yii ^ C»-(i) — 1 ^r(d— k— 1)- (5) 

Let t be some other term, not the diagonal term, in the determinant of Y. Thus, 

^ — Y d — k ,-K(d—k) ' ' ' Y dtir ( d ) ( 6 ) 

where tt is a permutation of {d — k , . . . , d}. Let j be the first index such that 7 r(j) ^ j, so 
that j lies between d — k and d. (Such a j must exist because t is not the diagonal term). 
Note that for this particular j , 7 r(j) > j, because if 7 r(j) < j then the permutation would 
repeat a column index from a previous factor. Thus, by the bounds in the previous section, 
we know that 

l^7. 7r (7)l — e Rj) “h e r(d-fc- 1)- (7) 

For j' > j, we get the crude upper bound 

|F)',7r(j')| < 1 + 2e 0 (8) 

because vy and uj-jt-i both lie in the rectangle [— e 0 , 1 + tof. 

Combining (4) and (5) yields 

0 > Yd-k4-k • • • Yj-ij-i ■ (er(j)-i — e r (d_jt_i)) • • • (e r (ci)-i - ^r(d-k-i))- 

Since the r’s are in monotone nonincreasing order, we know that r(d — k — 1) > r(j) > 
r(j + 1) > • • • > r(d). Thus, each factor in the preceding inequality from the jth onward is 
greater than or equal to — e r (j). Thus, the preceding inequality implies 

9 > Y d - k ,d-k ■ • ■ Yj-ij-i ■ ( e r(i)-i - e r (j)) d ~ J+1 . (9) 

On the other hand, if we substitute (7) and (8) into (6), we obtain: 

1*1 < U-M-* • ' ' Yi-U-i ■ (2e rW ) • (1 + (10) 
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Combining (9) and (10), we conclude that we can be assured that 


provided that 


*>|<|. ((* + 1 ) 1 - 1 ) 


( 11 ) 


2e r(j) • (1 + 2e 0 ) d 3 • ((* + 1)! - 1) < (e r (j)-i - tr(j)) d 3+1 (12) 

If (11) holds, then the sum of the other (k + 1)! — 1 terms of the determinant cannot cancel 
out 6 and the determinant of Y cannot be zero. This in turn implies that Bd-k - i cannot 
have a close point, i.e., I > d — k. I 

7 Choosing the Cj’s 

We now explain how to choose the e,-’s to satisfy (12). First, we want to get rid of j as a 
free index. Recall that j > d — k. Thus d — j < k. Therefore, if the following bound holds 
for each r = 0, . . . ,d — k — 1: 

2e r • (1 + 2co) fc • ((* + 1)! — 1) < (e r -i — Cr) fc+1 (13) 

this would imply (12). Furthermore, by making the inequality stricter, we can get a positive 
lower bound on det Y. For example, we could require more strongly that 

2e r • (1 + 2eo)* • (k + 1)! < (e r -i — Cr) fc+1 (14) 

which would give a lower bound of 6/(k + 1)! on det Y. 

To define the e r ’s, assume d and k are fixed. We now solve (14) with r = 0 to obtain Co 
(recalling that e_i is taken to be 1 — eo)- Next, we would find t\ from eo, and so on, up to 
£<*-*- 1- 

It does not seem possible to solve (14) in closed form to get an equation for e T in terms of 
e r _ i. We can get some crude estimates for general d, fc, r as follows. Let us assume first that 
e r < O.Se,—! for each r. Then the right-hand side of (14) has a lower bound of (0.5e r _i) fc+1 . 
Thus, (14) is satisfied by the recurrence: 


♦ 


J 


( 0 . 5 e f _ 1 )*' t ~ 1 
2 • (1 + 2e 0 ) fc •(* + !)!* 


To make further estimates, let us overestimate eo = 0.5 in the factor (1 + 2e 0 ) fc : 


( 0 . 5 e r - 1 ) fc+1 
2 • 2 fc • (* + 1)! 


Q2k+2 . (jfc+1)!* 
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To start this recursion off with r = 0, we can estimate e_i as 0.5 (recall that = 1 — Cq\ 
a suitable lower bound on e_i would be 0.5 since to < 0.5). Then the recurrence can be 
exactly solved starting with e_i = 0.5 to obtain 

tr _ q 5 (3fc+3)(/:+ir+(2fc+2)(fc+l) r - 1 +-+(2fc+2)(fc+l) 0 . + l)!](*+ 1 ) r +(*+ 1 ) r_1 +'"+(*+ 1 )°. 

The maximum value of r is d — k — 1 and the maximum value of k + 1 is d. Thus, we 
obtain as a worst-case estimate e r > 2~ 2 ° (<ll0 * d) . This bound is double-exponentially small in 
d because the second-level exponent is polynomial in d. 

Another approach to (13) is to solve numerically for particular values of d and k. We 
find a large value of e 0 satisfying (13) using the fact that e_i = 1 — eo- Then we find a large 
value of ci, and so on. We have carried out this computation for d = 3 and arrived at the 
following values that satisfy (13): 

For the case d = 3, k = 1, eo = 0.15, t\ = 0.007, €2 = 63 = 0. 

For the case d = 3, k = 2, eo = 0.05, d = £2 = £3 = 0. 

8 An aspect ratio bound 

Now, finally, we can combine the ingredients to obtain a bound on the aspect ratio of any 
simplex in our triangulation. According to Lemma 3, it suffices to obtain a positive lower 
bound on the determinant of the matrix 


W = 


Vi - v 0 
v 2 — v 0 


Vd — v 0 


This matrix can be simplified because we have now proved that Vo = (0, . . . ,0). This 
follows from Lemma 4: since d — k > 1, Bo does not have a close point. Thus, we need a 
lower bound on the determinant of 


t>i 



L J 


Following the analysis of Section 5, we conclude that W can be written in the form 


W = 


L 0 
X R 


where L is an (/ — 1) x (/ — 1) unit lower triangular matrix. Therefore, det(W) = det (R), so 
it suffices to obtain a positive lower bound on the determinant of the (d— l+l) x (d— / + 1) 
matrix R. 
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The worst case bound on the determinant of R comes when l = d — k, so that R is a 
(k+ 1) x (fc + 1) matrix. This matrix R has a very similar structure to the matrix Y analyzed 
in Section 6. In particular, the ( i,i ) diagonal entry of R is at least e r (,)_i, whereas the (i,j) 
superdiagonal entry is at most e r ( t ) in magnitude. The subdiagonal entries are at most 1 + e 0 - 
These bounds are in fact better than the bounds on Y, in other words, the diagonal entries 
of R have better lower bounds than the corresponding diagonal entries of Y, and similarly, 
the off-diagonal entries have better upper bounds. Thus, the determinant of R is at least 
l/(k + 1)! multiplied by the diagonal term, because the same result held for Y if (14) holds. 
The worst case is when all the diagonal entries of R are td-k- i- Thus, a lower bound on the 
determinant of W is + 1)! Recall that td-k-i is double-exponentially small in d. 

Therefore, raising this to the k + 1 power and dividing by (fc + 1)! yields a lower bound on 
the determinant that is also double-exponential in d. 


9 Open questions 


Here are some open question raised by this work: 


1. Can the bounds be improved? Is the aspect ratio bound inherently double exponential 
in d (i.e. are there lower bounds)? How much can the actual numerical bounds be 
improved for d = 3? 

2. Can this work be generalized to the case that F is a curved manifold? This would be 
a first step towards generalizing the QMG algorithm to handle curved surfaces. In the 
curved case, it is presumably necessary that the grid of boxes be sufficiently small that 
F bends only a small amount in each box. 


3. Is it possible to triangulate a grid of boxes cut by a hyperplane in such a way that all 
the dihedral angles of all simplices are bounded by 7t/2? Dihedral bounds of tt/ 2 are 
important for some applications (see e.g. [13]) but it is not known how to triangulate 
even the special case of convex three-dimensional polyhedra and obtain such an angle 
bound. (The latter problem was solved in two dimensions by [1]; see also [2], [4] and 
[ 8 ].) 
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